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ABSTRACT 






A comprehensive numerical model has been developed to investigate 
the characteristics of an impu I s i ve I y-started , asymptotica I ly-steady, 
uniform flow about a circular cylinder. The model is based on the 
red i scret i zation of the shear layers, wake-boundary- I ayer interaction, 
and the dissipation of vorticity by turbulence. 

The forces acting on the cylinder, rate of vorticity flux, 

Strouha I number, oscillations of the stagnation and separation points, 
longitudinal and transverse spacing of the vortices, and the base 
pressure have been calculated and shown to be in conformity with those 
obtained experimentally. 

The model has been used to predict the character i st i cs of hydro- 
elastic oscillations of a cylinder in the range of synchronization. 

The numerical experiments shed considerable light on the interaction 

between the fluid motion in the wake and the dynamics of the body. 

* 

An extensive sensitivity analysis has been carried out to determine 
the stability of a I I the parameters and hence the stability of the 
numerical model itself. 
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I . INTRODUCTION 



A. FLOW SEPARATION 

Real fluids flowing over bluff bodies separate under the influence of 
adverse pressure gradient as momentum is consumed by wall shear. Separa- 
tion is best understood in terms of its consequences. In fact, it was the 
keen observations of separation in diffusers that led Prandtl to introduce 
the concept of transition or boundary layer theory. Even though the 
boundary layer theory has revolutionized fluid dynamics, the phenomenon 
which gave impetus to its discovery still remains an enigma. 

Through the years various character i st i cs of separated flow about 
bluff bodies in general and about circular cylinders in particular have 
been partially understood through observations, measurements, and analyses. 

The separat i on point may be mob i I e as on a ci rcu I ar cy I i nder or f i xed 
along a line as on a wedge or blunt-based body. When a separation point 
is mobile, the pressure decreases up to the point of maximum velocity 
(slightly upstream of separation) and then increases up to the point of 
separation. When the separation point is fixed, the pressure continues 
to decrease up to the point of separation. Consequently, there are signi- 
ficant differences in the separation of flow about a circular cylinder 
and a sharp-edged body. In either case the velocity d i str i bution, rate 
of vorticity flux, and the amplitude and frequency of oscillation of the 
separation points are not independent of the wake in the downstream side 
of the body. It is this interaction that makes the calculation of sepa- 
rated flows exceedingly complex. In other words, one needs the evolution 
of the boundary layer on the forebody (upstream of separation) of the 
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cylinder in order to determine the instantaneous position of the separa- 
tion points. But, the characteristics of the boundary layer are, in turn, 
strongly affected by the fluid motion in the wake. Evidently, one needs 
some kind of an interactive mechanism whereby the effect of the wake on 
the boundary layer, and hence separation, can be given proper considera- 
tion. In mathematical terms one can state that the unseparated potential 
flow, or the d'Alembert flow, does not even constitute a first approxima- 
tion to the separated real flow. Partly the realization of this fact, 
partly the difficulties associated with the integration of the Navier- 
Stokes equations across a singularity (separation point), and partly lack 
of sufficient computational capability have prevented the prediction of 
the kinematic and dynamic characteristics of separated flows about bluff 
bodies. In the following, the basic characteristics of flow about circu- 
lar oylinders are described with special emphasis on those aspects which 
prepare the reader for the understanding of the numerical model developed 
during the course of this work and the assessment of the similarities and 
dissimilarities between the predictions and observations and measurements. 

B. VORTEX SHEDDING 

The separation of flow gives rise to alternate vortex shedding. This 

phenomenon is intrinsic to the flow itself and has been known at least 

since the times of Leonardo de Vinci. In 1878, almost exactly 100 years 

* 

ago, Strouha I discovered that there is a relationship between the frequency 
of vortex shedding, the velocity of flow, and the diameter of the cylinder. 
This relationship, which is now known as the Strouhal Number, St, is 
given by St = fD/U. Strouhal 's experiments gave an average value of St = 
0.185. During the course of numerous experiments since then, it has been 
experimentally demonstrated that: (i) the Strouhal number is dependent on 



16 



the Reynolds number and, as first suggested by Lord Rayleigh, must be 
written as f = St(Re)U/D; (ii) the vortex shedding frequency is only an 
average quantity for a given flow and that there are always secondary and 
tertiary frequencies present; and finally, (iii) there are regions of flow 
such as the region of drag crisis where it is quite difficult to identify 
a particular frequency. It is evident from these facts that in some re- 
gions of flow one must speak of the spectral content of the wake fluctua- 
tions in assessing the value of the Strouhal number and the forces exerted 
on the cylinder. The foregoing discussion of the Strouhal number raises 
two questions: firstly, is it possible to devise a universal Strouhal num- 
ber which would be applicable to a I I two and three-dimensional bluff 
bodies, and secondly, is it possible to theoretically predict the Strouhal 
number through, for example, the perturbation of the d'Alembert flow. An 
attempt along these lines has been made by Sacksteder C I H who found St = 
0.2028 for an infinitely large Reynolds number. The reasoning behind 
this prediction and its validity will not be discussed here further. 

The former question of the possibility of devising a universal Strouhal 
number has been of interest to various researchers. Table I, shown below, 
lists those definitions together with their originators. 



TABLE I 



St = 


Va 1 ue 


Body Type 


Or i g i nator 


fyh'/UK 


0.163 


two-d imens Iona I bod i es 


Roshko [2] 


f % /u 

V V 


0.181 


two-d i mens i ona 1 bod i es 


Bearman 


W K " J 


0.190 


ax i symmetric bodies 


Ca 1 vert C4j 


f vV K ’ U 


0.163 


two-d imens iona 1 bod i es 


Simmons L5j 



In Table I, h’ represents the wake width as determined from Roshko’s 
notched- hodograph theory C23; & v , lateral spacing of the vortices as 
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determined by Kronauner f s minimum drag criteria as reported by Bearman C33 

i , the wake width at a distance of the wake establishment region; and K = 
w 

/f-C * ~ where C . is the base pressure coefficient. All these definitions 

pb pb r 

require either the measurement or the solution of one or more characteris- 
tics of the wake. Thus, their capability to predict a vortex shedding fre 
quency is quite limited. 

It is advantageous to note at an early stage those factors which in- 
fluence the Strouhal number. Firstly, small amplitude transverse oscilla- 
tion of cylinders or cables distinctly regularize the vortex shedding and 
hence the Strouhal number. For example, the experiments conducted at the 
National Physical Laboratory CSU have shown that the Strouhal number even 
in the drag crisis region remains perfectly constant and equal to about 
0.21. Secondly, the vortices do not peel off from a cylinder as perfectly 
two-dimensional line vortices. As it was first shown by Humphreys 
there is a cell pattern separated by a correlation length. The longer the 
cylinder the larger the number of cells. The consequences of the forma- 
tion of such cells are that there are phase shifts in pressure along the 
cylinder and hence a reduction in the integrated transverse force. In 
other words, the transverse force calculated through the integration of 
the pressure distribution at a given section is not identical to the aver- 
age transverse force determined experimentally. All numerical analyses 
which model the complex three-dimensional flow as a two-dimensional flow 
must necessarily predict a transverse force or lift coefficient somewhat 
larger than those found by measurements. Other consequences of the cell 
formation concern the strength of the vortices. The vortices are not 
necessarily oriented in a plane normal to the flow and the measurements 
(flow visualization or electro-optical techniques) yield only that 
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component of vorticity which is normal to the plane of flow. Additional- 
ly, experimental results necessarily depend on the I ength-to-d iameter 
ratio, end conditions, intensity and scale of turbulence in the ambient 
flow, the rigidity of the cylinder mounting, and of course on Reynolds 
number and the surface conditions of the cylinder. Most recent experi- 
ments (Shih and Hove C8H) have shown that roughness affects the base 
pressure in the same manner as the small transverse oscillations (beyond 
a threshold limit of about 0.10 diameters). Uniform roughness yields a 
perfectly uniform base pressure and removes the cell pattern. 

C. RESISTANCE IN STEADY FLOWS 
I . In-Line Forces 

The understanding of the fluid forces exerted on bluff bodies 
comes primarily from measurements. Of all bluff bodies, the circular cy- 
linder has been subjected to the largest number of investigations. Those 
experiments which constitute milestones in the determination of the magni- 
tude and characteristics of resistance are the works by W iesel sberger C^H, 

Fage and Johansen ClOj, and Rhosko [l l], A plot of the drag coefficient, 

2 

0^ defined by = 2F/pLDU , is shown in Fig. I as a function of the 

Reynolds number defined by Re = UD/v. It is a well known fact that for 

4 5 

Reynolds numbers larger than about 10 and smaller than about 2x10 , the 
drag coefficient for a circular cylinder remains essentially constant at 

5 

a value of about 1.2. For Re > 2x10 , the shear layers emanating from 
the separation points become turbulent. Since a turbulent flow can ex- 
change larger momentum between the inner and outer flow ( i . e . , between the 
boundary I ayer and externa I flow), it can sustain a larger adverse pressure 
gradient and even reattach to the cylinder. The reattachment of the tur- 
bulent shear layers gives rise to the formation of so called separation 
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Drag coefficient versus Reynolds number for a smooth circular cylinder. 



bubbles and to the displacement of the attached flow further downstream 
from the first separation point. The reattached flow cannot continue to 
remain attached indefinitely and reseparates at an angle of about 130 de- 
grees from the front stagnation point. Thus, the transition of the flow 
in the shear layers from a laminar to a turbulent state results in a con- 
siderably smaller wake and hence in a considerably reduced drag coefficient. 
This particular flow regime is commonly called the drag crisis region or 
the critical region. As the Reynolds number increases, parts of the bound- 
ary layer adjacent to the first separation point become turbulent and the 
second separation point begins to move gradually upstream. This region is 
called the transcr i t i ca I region. At higher Reynolds numbers, the entire 
boundary layer becomes turbulent and the drag coefficient reaches a nearly 
constant value of about 0.6. There is considerable scatter in the data 
in the supercr i t ica I region because of the difficulties associated with 
flows at such high speeds. 

The present knowledge of the drag coefficient for a smooth cylin- 
der is limited to Reynolds numbers smaller than about 10^. There is great 
practical interest in the determination of force coefficients, i.e., drag 
coefficient, lift coefficient, and the Strouhal number for projects such 
as ocean thermal energy conversion (OTEC), large cooling towers, high rise 

structures, etc. A careful cons i derat ion of the various possible methods 

8 1 0 

to obtain data at Reynolds numbers in the order of 10 or 10 shows that 
there are practically insurmountable difficulties. The existing wind and 
water tunnels are not capable of yielding such high Reynolds numbers with- 
out the adverse effects of blockage and compressibil ity. The pop-up and 
drop-down experiments with large cylinders are constrained by cavitation, 
weight, and stability conditions. It might be possible in the future to 
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close the high Reynolds number gap by building large cryogenic wind tun- 
nels operating with freon or liquid helium. 

The drag coefficient depicted in Fig. I may be influenced by a 
number of factors such as roughness and turbulence. Numerous investiga- 
tors, notably Fage and Warsap [12], Achenbach D3H, and Guven El4U, have 
clearly demonstrated that the larger the surface roughness, the sooner is 
the occurrence of the critical region and the higher the dip in the drag 
coefficient (see Fig. 2). Turbulence plays a similar role in the evolu- 
tion of the drag coefficients. It should be noted from Fig. 2 that the 
drag coefficient for a rough cylinder at supercr i t i ca I Reynolds numbers is 
about twice that for a smooth cylinder. Furthermore, the experiments by 
Guven L 1 4U and Szechenyi 0 5] have shown that at sufficiently high 
Reynolds numbers there is a Reynolds number independent region. This ex- 
perimental fact has been termed the independence principle and used (Shih 
and Hove CO]) to extrapolate the drag coefficient data to Reynolds numbers 

g 

in the order of 10 . If the validity of the independence principle can 
be demonstrated unambiguously, it may obviate the need for difficult ex- 
periments in cryogenic tunnels. 

2. Transverse Forces 

Lift or transverse-force coefficient data are far less numerous and 

exhibit a great deal of scatter. Figure 3, which represents most of the 

existing lift coefficient data, shows that the lift coefficient, defined 

2 

by = (amplitude of lift force)/0.5pLDU , may vary from 0.2 to 1.0 in the 
subcritical Reynolds number range. In the critical range it is impossible 
to define a lift coefficient without resorting 1o its root-mean-square (rms) 
values. It is for this reason that some of the data appearing in Fig. 3 
are in terms of the rms values of the lift coefficients. For a more 
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Drag coefficient versus Reynolds number for a circular cylinder with distributed roughness Cl4[]. 
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Experimentally obtained values of lift coefficient versus Reynolds 
number for a circular cylinder C I fc>3 • 



detailed discussion of the spectral content of lift in the critical region 
see Bub I i tz 0 7], 

The reasons associated with the large variations of the lift co- 
efficient are directly related to the mobility of the separation points, 
formation of vortex cells, lack of spanwise correlation, end effects, tur- 
bulence present in the ambient flow, etc. Since each experiment is carried 
out with a different I ength-to-d iameter ratio, blockage ratio, and end 
conditions, the number of vortex cells over the cylinder and the effect of 
the flow between the end of the cylinder and the tunnel wall give rise to 
widely different lift coefficients. 

The effect of roughness on the lift coefficient is to reduce its 
scatter and to increase its magnitude. As noted earlier, roughness dis- 
tinctly regularizes the vortex formation and the base pressure thereby 
eliminating the effect of phase shifts which cause the said scatter in 
lift. Data supporting the above experimental facts will not be reported 
here for sake of brevity. For additional details the reader is referred 
to the work by Szechenyi Cl 5], 

D. RESISTANCE IN TIME-DEPENDENT FLOWS 

There are an infinite number of time dependent flows and it is impos- 
sible to discuss their character i st i cs in any unifying manner. One might 
obtain a clearer insight into such flows by considering only those cases 
which are relatively more manageable for numerical or laboratory experiments. 
Among such flows which have been subjected to extensive numerical analysis 
and experiments are the impulsively started flows, uniformly accelerating 
■flows, and harmonically oscillating flows. Obviously, no distinction is 
made here between the said fluid flows about a body at rest and a body 
undergoing similar time dependent motions in a fluid otherwise at rest. 
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I . Impulsively Started Flows 



The impulsively started flow has been by far the one most investi- 
gated numerically through the use of the Nav i er-Stokes equations. These 
investigations (Collins and Dennis and Yang [I9j) have been confined 

to the very early stages of the motion during which the separation points 
move from the rear stagnation point to about 109 degrees. An exception 
to the foregoing is the work reported by Thoman and Szewczyk C20X They 
have solved the Nav i er-Stokes equations through the use of a suitable 
finite difference technique for various Reynolds numbers and found that 
the drag coefficient overshoots for a relative flow displacement (normalized 
time) of Ut/c = 5 and then gradually reaches steady state values. 

Only two separate investigations of the impulsively started flow 
have been reported. The first is that carried out by Schwabe C2 I J upon 
Prandtl T s suggestion and the other by Sarpkaya [ ' 22 , 23j. Schwabe ! s experi- 
ments were conducted with a 9 cm diameter cylinder at a Reynolds number 
of about 600. His drag coefficient, calculated indirectly through the use 
of motion pictures, rapidly reaches a value of about 2.0 and remains con- 
stant thereafter. Schwabe's findings are not in conformity with either 
the numerical calculation or with the experimental results obtained by 
Sarpkaya. Figure 4 shows three sets of representative data at Reynolds 
numbers from 16000 to 35000 for the drag coefficient. Clearly, the very 
initial instants of the motion cannot be analyzed because of the problems 
associated with the generation of an impulsively started flow. At Ut/c 
slightly greater than zero, the force acting on the cylinder is almost 
entirely inertial. As the separation sets in at the rear stagnation point 
and rapidly moves upstream, the drag also increases rapidly and reaches 
a value of about 1.5 at Ut/c between 4.0 and 5.0. Then C D drops rapidly 
to the steady state value for smooth cylinders in steady uniform flows at 
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Fig. 4. Drag coefficient versus normalized time for impulsively-started flow about 
a circular cylinder C23], 



Ut/c larger than about 7.0. Detailed flow visualization experiments by 
Sarpkaya [22 H have shown that the reason for the overshoot i n Cp is the 
symmetric development of two vortices in the cylinder wake and the rapid 
accumulation of vorticity in each vortex. It must be noted that for such 
a symmetric vortex configuration, the only means by which the vorticity 
may be cancelled are the exchange of vorticity between the two vortices 
and the cancellation of vorticity with oppositely-signed vorticity, so 
called counter-vort i c i ty , generated in the rear boundary layers. In other 
words, the rate of accumulation of vorticity in the vortices is consider- 
ably larger than that at any other time. Other vorticity cancellation 
mechanisms which come into play during the separation of the vortex from 
its feeding sheet, as proposed by Gerrard [24j, do not play any role dur- 
ing the growth of the symmetric vortices. Thus, the rapid accumulation 
of vorticity gives rise to the drag overshoot. As the asymmetry sets in, 
vortices shed a I ternat i ng I y and the flow gradually becomes asymptotically 
steady. The possible mechanisms which cause the asymmetry in an otherwise 
symmetrical flow about a symmetrical body are discussed in greater detail 
later. Suffice it to note that one of the fundamental facts of fluid 
flows is that symmetric causes do not give rise to symmetric effects, as 
most eloquently philosophized by Birkhoff [25^. 

Following the onset of asymmetry, the cylinder experiences a 
transverse force, as shown in Fig. 5, in terms of the lift coefficient and 
the normalized time Ut/c. The lift force fluctuates with the Strouhal fre- 
quency. The shedding of each vortex gives rise to a small fluctuation in 
0^. The frequency of this fluctuation is twice the Strouhal frequency 
since two vortices (one from top and from bottom) are shed during a given 
cycle of I i ft osci I I at ions. 
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Lift coefficient versus normalized time for impulsively-started flow about 
circular cylinder C23H. 



2. Other Ti me-Dependent Flows 



As noted earlier, the uniformly accelerating flows and harmonical- 
ly oscillating flows have been subjected to numerous analytical and experi- 
menta I investigations. Among those notable are the investigations carried 
out by Iverson and Balent [ 262 , Keim [. 272 , Hamilton [ 282 , Odar and Hamilton 
[ 292 , Sarpkaya and Garri son. [30], and Sarpkaya C3 I H- The basic question 
has been and continues to be the formulation of a resistance equation which 
takes into consideration the effects of acceleration, velocity, and the 
history of the motion. All investigations have shown that resistance in 
time-dependent flows in general and in harmonically oscillating flows in 
particular cannot be considered as a juxtaposition of resistance in steady 
flows at correspond i ng instantaneous Reynolds numbers. Stokes, in his 
celebrated paper on the sinusoidal oscillation of a sphere in a fluid 
otherwise at rest, has shown that the fluid resistance is comprised of a 
linear velocity-dependent force and an acceleration-dependent inertial 
force. The coefficients for both forces depend on the Reynolds number and 
the frequency parameter D /vT where T is the period of oscillation. Ob- 
viously, Stokes T analysis is valid for only very small, as yet unspecified, 
Reynolds numbers. Inspired partly by Stokes* contribution and partly by 
practical cons i derat ions, Morison [522 proposed a resistance equation which 
may be written as 

F = I pA p C D u l u l + pV C M 3T n> 

in which U represents the instantaneous velocity; A , the projected area; 

V, the displaced volume; and and C^, the drag and inertia coefficients 
respectively. For harmonically oscillating flow about a stationary cylin- 
der where the fluid velocity is IJ = -U^ cos 0 (here 0 = 2frt/T) Eq. (I) 
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reduces to 



7 p0L U M 2 



- C D |cos 0 1 cos 6 + ^ T/D C M Sin 0 



( 2 ) 



The essence of this equation is that the resistance in harmonic flow is 
assumed to be equal to a linear combination of a ve loc i ty-squared dependent 
drag force and an acceleration-dependent inertia force. In the midst of 
i nsuff ic ient I y clear insight into the kinematics of the complex flow, it 
has not yet been possible to demonstrate the reasons for the differences 
between the measured force and that calculated from Mori son's equation for 
U^T/D values in the vicinity of 15. Additional experiments and numerical 
analysis, possibly through the use of the discrete vortex model, might shed 
further light on the subject and into the overa I I understanding of resis- 
tance in time-dependent flows. It is evident from the discussion of both 
the steady and unsteady flows that analytical or numerical methods are 
urgently needed with which relatively inexpensive investigations can be 
conducted to clarify the complex wake-boundary- I ayer interaction mechanism. 
Of all the possible flows which might be considered as a candidate for such 
an analysis, the impulsively started flow is certainly the one most manage- 
able and most fundamental to the understand i ng of ail other flows. 



E. METHODS OF NUMERICAL ANALYSIS 

With the advent of high speed computers, numerical methods began to 
play a greater role in fluid dynamics. Currently four major methods of 
analysis are in use: finite difference, finite element, marker-cell tech- 

nique, and the discrete vortex method. 

Finite difference methods have been used for the solution of flow 
about a ci rcu I ar cy I inder for relatively small Reynolds numbers (Re - 100) 
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either through the use of the Nav i er-Stokes equations as they are commonly 
written or through the more frequently used vorticity equations. Ingham 
[33] calculated the steady flow past a circular cylinder at Re = 100. 

Thoman and Szewczyk [20], as noted earlier, carried out a similar numeri- 
cal analysis of flow over stationary and rotating cylinders at Reynolds 
numbers as large as 10 . There are a number of such examples in the 
I iterature but they are not cited here because of their lack of relevance 
to the subject under consideration. 

The finite element method, which came into prominence first in the 
area of structura I analysis and then permeated through heat transfer and 
fluid mechanics, has been primarily used for unseparated flows about 
bluff bodies. As far as fluid mechanics is concerned, the outstanding ex- 
amples of such applications may be found in two recent volumes of work 
dealing with finite elements in fluid mechanics [34]. • Bretanow and Ecer 

[35] made an attempt to calculate the separated flow about a circular 
cylinder at small Reynolds numbers but the results have been i nconc I us i ve. 
Applications of the finite element method to flows about foils where the 
flow remains attached have yielded results in conformity with experiments 
(Bretanow and Ecer [35]). 

The marker-cell technique developed by the researchers at the Los 
Alamos Scientific Laboratory during the past 20 years uses the finite 
difference form of the Nav i er-Stokes equations and then adjusts the pres- 
sure in each cell to insure that the equation of continuity is satisfied. 
Numerous ingenious techniques have been developed by Fromm and Harlow 

[36] and Nichols and Hirt [37] for the efficient adjustment of pressure 

in each cell. Such techniques will not be dealt with here. The researchers 
of Los Alamos Scientific Laboratory have applied the marker-cell technique 
to the flow past cylinders, spheres, rectangular blocks and to various other 
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free-surface flows. These applications dealt with very low Reynolds num- 
bers (in the order of 100). Even though the technique has considerably 
increased the understanding of the kinematics of the wake, the precise 
nature of the time-dependent force acting on any one of the bodies cited 
above remains obscure. 

The discrete vortex model, which is included as one of the four major 
models, is the subject of extensive discussion in the next chapter. There- 
fore, the remainder of this section is confined to a discussion of the suc- 
cesses and shortcomings of the first three methods. 

All numerical models are by their very nature approximate. They satis- 
fy certain critical conditions only at discrete time steps. Secondly, they 
requ i re mesh or element sizes which vary from one region of flow to another 
to insure that the character istic gradients of the flow (such as velocity 
or pressure) do not significantly or abruptly change over any one element. 
For example, the finite difference method when applied to the viscous flow 
over a circular cylinder requires that the size of the elements on and near 
the cylinder be considerably smaller than the boundary layer thickness. 
Otherwise, the evolution of the boundary layer, the positions of the sepa- 
ration points, and the evolution of the wake cannot be accurately calcu- 
lated. Since the boundary layer thickness for laminar flow is proportional 
- 1 /2 

to Re , one realizes that only for relatively small Reynolds numbers can 
a sufficiently large boundary layer thickness and a correspond i ng I y large 
mesh size be obtained. This fundamental difficulty has prevented the 
finite difference, finite element, and the marker-cel I technique from making 
significant progress into the calculation of steady or time dependent flows 
about bluff bodies at Reynolds numbers of practical interest. None of the 
three techniques cited above has been applied to the solution of flow prob- 
lems associated with oscillating bodies. The reasons for this are rather 



33 



self-evident since the solution requires not only a moving grid but also 
a grid whose mesh size must be varied with time to minimize the effect of 
sharp velocity and pressure gradients. Finally, it should be noted that 
all three numerical techniques present "stiffer" solutions and introduce 
an artificial numerical viscosity whose influence cannot be readily 
assessed. Not to be entirely critical of the said methods, for there are 
corresponding shortcomings of the discrete vortex model, one must note 
that in the wide spectrum of Reynolds numbers there will always be some 
interest in low Reynolds number flows. Thus, such numerical techniques 
will continue to play some role in separated flows where the discrete 
vortex model is not advantageously suited. It therefore appears that no 
one model will stand out as superior and un i versa lly applicable and that 
all methods will continue to complement each other. 
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II. FUNDAMENTALS OF THE DISCRETE VORTEX MODEL 



A. INTRODUCTION 

The discrete vortex model (DVM) is a potential flow representation of 
the shear layers emanating from the separation points. Any separated vis- 
cous flow is completely specified if the spatial and temporal distribution 
of vorticity in the entire flow field is known. For reasons described in 
the previous chapter, it has not yet been possible to achieve this objec- 
tive through any exact or approximate methods. It has been observed and 
experimentally verified by Fage and Johansen 0 Oj that the shear layers 
at sufficiently high Reynolds numbers are quite thin and the vorticity is 
confined into regions of tightly spiralled vortex sheets. This fact alone 
invites one to divide the shear layer into a number of small segments and 
concentrate the vorticity of each segment into a line vortex situated at 
the geometric center of the segment. This idealization permits one to 
use the powerful complex function theory to determine the kinematics and 
dynamics of the flow. Evidently, the validity of the results can only 
be judged in light of experimental facts. 

In the following, only the gross features of the method are described 
and the details are presented following an historical resume. 

At a given time the separated shear layers are represented by vortices 
situated at the center of each segment. The potential flow theory permits 
one to calculate the velocity of each vortex induced at its center by all 
other vortices and the ambient flow. Then the vortices are convected for 
a short time interval At through the use of a suitable convection scheme 
to new positions. It is then necessary to introduce new vorticity into 
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the flow at or near the separation points to account for the vorticity 
generated in the boundary layer. The introduction of new vortices (nascent 
yortices) represents an addition to each set of vortex arrays. Then the 
process is repeated. 

In spite of its apparent simplicity, there are a number of problems 
which must be resolved in carrying out the said calculations: (i) the vor- 
ticity introduced at each time step must be related to the character i st i cs 
of the boundary layers upstream of the separation points; (ii) one must 
insure that the point vortices continue to represent the vorticity at each 
segment in accordance with the fundamentals of hydrodynamics; (iii) allow- 
ance must be made for the interaction of oppositely-signed vorticity and 
the decay of vortices; (iv) methods must be devised to limit the number of 
the point vortices to minimize the computation time without affecting the 
kinematics and dynamics of the flow; and (v) procedures must be established 
for the flow to retain its inherent flexibility without introducing 
artificial stiffening effects (e.g., to allow for the mobility of stagna- 
tion and separation points). 

The foregoing conveys the idea that the discrete vortex model is and 
will continue to rely on a strong interaction between potential and vis- 
cous flow theories. Neither is capable of providing a solution alone but 
together they can make inroads towards describing the basic features of 
the separated flow provided that sufficient ingenuity is incorporated into 
the said interaction. Lastly, it must be noted that all the numerical 
methods described above depend on the availability of a high-speed digital 
computer. 

B. HISTORICAL DEVELOPMENT 

The DVM has been in existence for the past 46 years since its incep- 
tion by Rosenhead [38j in 1932. Rosenhead investigated the so-called 
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"Helmholtz instability" of a surface of discontinuity, separating two 
fluids of equal speed and density but moving in opposite directions. His 
primary objective was to give the shear layer between the two fluids an 
initial disturbance by situating the shear layer along the curve y = 0. I 
sin (2 ttxA) and then observe the evolution of the sheet with time. To ac- 
complish this, Rosenhead represented the shear layer by point vortices 
placed along the sheet at intervals such that their projections along the 
x-axis were equ i-d i stant. Each vortex was assigned a strength of UX/6n 
where n is the number of vortices per wave I ength, X, and U the' speed of 
the two fluids. Keeping the vortex strengths constant, Rosenhead followed 
the deformation of the sheet by convecting the vortices with a simple 
Eulerian integration of = rj + q^At where r ^ is the position vector 

of the j-th vortex; q J , the computed velocity vector; and At, the time 
interval. Using n = 2, 4, 8, and 12, Rosenhead observed that in each case 
the sheet rolled up smoothly into concentrated clusters of vortices posi- 
tioned at intervals of X along the x-axis. Rosenhead f s results, which were 
computed using a desktop calculator, are shown in Fig. 6. The tendency of 
a distribution of point vortices to roll up in a manner consistent with 
observations is clearly evident. 

Rosenhead r s results remained undisputed until 1959 when the computa- 
tional power of the digital computer was applied to the analysis. Birkhoff 
and Fisher [ ’39 U repeated Rosenhead’s analysis using a Runge-Kutta integra- 
tion scheme, 22 vortices per wavelength, and smaller time intervals. In- 
stead of improved results they observed that the vortices took irregular 
paths and no longer rolled up smoothly. Birkhoff and Fisher concluded that, 
given sufficient time, an originally organized array of equ i -strength point 
vortices would eventually become random and that the validity of represent- 
ing a shear layer by point vortices is questionable. 
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Fig. 6. Example of Rosenhead's 
calculations C383. 
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Shortly thereafter, Hama and Burke C40H undertook a similar computer 
study and verified the observations of Birkhoff and Fisher. However, Hama 
and Burke noted a conceptual error in Rosenhead's original method of plac- 
ing the discrete vortices along the shear layer. Specifically, if the 
initial vortex sheet was to have a constant vorticity per unit distance 
along its length and be represented by equ i-strength vortices, then the 
vortices should be placed at equal distances measured along the sheet vice 
at equal x increments. Using this arrangement Hama and Burke repeated 
their calculations and found -that the onset of irregular vortex motion was 
significantly delayed. Eventually, however, a random pattern did result. 
Thus, at this point it was hypothesized that the onset of instability 
could be delayed although the questions of why it occurred and how it 
could be control led were still a mystery. 

Without fully realizing the consequences of sheet instability (the 
result of which is discussed later), the DVM continued in its evolution 
relying on the assumption that meaningful results could be obtained up to 
the onset of irregular vortex motions. In 1962 Abernathy and Kronauer C4I]] 
applied the method to investigate the motion of two shear layers of oppo- 
sitely-signed vorticity separated by a distance h y and each given an 
initial perturbation as in Rosenhead's work. They observed that the vor- 
tices accumulated into clusters resembling the pattern of a vortex street. 
An example of their results is shown in Fig. 7. 

Encouraged by the demonstrated ability of the DVM to represent a vor- 
tex street, the next logical extension of the method was to analyze the 
flow past a bluff body. However, the inclusion of a solid body in the 
flow field necessitated new considerations. Specifically, the previous 
analyses had begun with a discretized sheet having a preassigned initial 
position and vorticity distribution along its length (usually taken as 
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Fig. 7. 



Example of Abernathy 

and Kronauer's calculations 

( h/X=0 . 34 ) [41]. 
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constant). The application of the DVM to flow past a bluff body required 
that the discretized shear layers be generated continuously by adding 
additional point vortices into the flow at regular time intervals and at 
locations representative of where the shear layers actually emanate from 
the body. The requirement for some method to account for vorticity genera- 
tion and introduction in an otherwise potential flow model posed a signi- 
ficant problem. Compounding the difficulty was the time dependent nature 
of both the generation rate and its point of introduction. 

The first DVM analysis to include a bluff body was performed by 
Gerrard C42J for a circular cylinder. Although a circular cylinder was of 
great practical interest it was also the most difficult one to analyze 
because of the unsteady vorticity generation and mobile separation points. 
Gerrard was forced to adopt a rather arbitrary method of introducing new 
vortices which severely influenced the results of his analysis, although 
it still represented a noteworthy application of the DVM. Noting the dif- 
ficulties associated with the mobile separation points and the vorticity 
introduction, Sarpkaya £43] introduced a mobile separation region repre- 
sented by an infinite number of vortex sheets of vanishingly small vorticity 
connected to the nascent vortex. This analysis was restricted to impul- 
sively started symmetric flow about a circular cylinder. In the ensuing 
years many investigators bypassed the difficulty introduced by the mobility 
of the separation points and dealt with bodies with sharp edges (e.g., 
inclined plates, blunt-based bodies, wedges, etc.) where the separation 
points are fixed. Despite the growing number of applications of the DVM, 
a number of disturbing difficulties remained unresolved and a full inter- 
action between the boundary layer, separation points, and the wake had 
not been established. These difficulties required the use of a number of 
nondisposable parameters. Consequently, the proper balance between 
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heuristic reasoning and computational formulation had to evolve through 
the recognition of the reasons underlying the basic difficulties. 

Briefly, the most important of these difficulties stemmed from the 
proximity of the vortices to each other and to the solid boundary, and 
from the necessity of representing the continuous process of vorticity 
generation at the singular separation points by a discrete process. The 
occasional proximity of two point vortices resulted in large mutually- 
induced velocities and hence the shooting of the two vortices rapidly 
away from the flow field. This was particularly true when the vortices 
ceased to be located at the geometric center of the sheet segment which 
they were supposed to represent. Furthermore, the spiralling of the sheets 
brought vortices on two adjacent spirals close to each other and gave rise 
to a traveling and growing instability. Thus, it was not possible to 
correctly represent the core of tightly spiralled sheets. The proximity 
of the vortices to the boundary and hence to their images inside the 
body resulted not only in similar instabilities but also in unrealistic 
local pressures. This, in turn, resulted in unnatural fluctuations in 
drag and lift forces. 

Faced with the problems just cited and the desire to enhance the 
applicability of the method, various i nvest i gators devised several numeri- 
cal cures. Chorin and Bernard C44] introduced a cut-off distance from 
the center of each vortex beyond which the vortices were to behave like 
potential vortices. For radial distances smaller than the cut-off distance 
the vortices were assumed to exhibit a solid body rotation. Their calcu- 
lations did not, however, result in an improvement of the analysis of 
separated flow about a circular cylinder. Moore E45j, in dealing with 
two vortex sheets emanating from a plane lifting surface, introduced the 
method of amalgamation. When a vortex sheet was spiralled tightly enough. 
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